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Abstract 



We investigate the neutrino oscillations of z/ M beam at the K2K experiment in 
the four-neutrino model with three active and one sterile neutrinos, and compare 
them with the oscillations in the three-neutrino model. In the four-neutrino case, 
the effect of the AmL SND scale of mass-squared difference, derived from the LSND 
experiments, occurs in the survival probability Piv^ — > Vy) in the range of Am 2 < 
0.004 eV 2 , where Am 2 is the relevant one to the K2K experiment and corresponds 
to the atmospheric neutrino mass scale. Once the probability Piy^ — > v^) is 
measured at the K2K, the allowed region of Am 2 would turn out to be broader 
in the four-neutrino model than the one in the three-neutrino model. 
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It has turned out that neutrinos have a certain amount of mass through the 
observations of the atmospheric neutrino anomaly |I] M. The anomaly can nat- 
urally be explained by the neutrino oscillation Q, along with the analyses of the 
solar neutrino deficit [Q. The oscillation, however, gives only the mass-squared 
difference among various species of neutrinos. 

At present, if the LSND experiment is included, three kinds of mass-squared 
differences are derived: Am 2 olar = (10~ n — 10~ 5 )eV 2 from the solar neutrino 
deficit with a large range of Am 2 , depending on the four solutions of the vacuum 
oscillation, and the Mikheyev-Smirnov-Wolfenstein (MSW) solutions in the matter 
with small- and large-angle mixings and the LOW one with relatively low mass- 
squared difference |)|], Am 2 tm = (1.5 — 5) x 10~ 3 eV 2 from the atmospheric neutrino 
anomaly with a large mixing angle of sin 2 2# atm > 0.82 interpreted as the — > u T 
oscillation, and Ai?Il SND = (0.2 — 2)eV 2 from the LSND experiments on — > v e 
and — > v e oscillations which is the only one positive evidence from the 
terrestrial oscillation experiments using the accelerators and reactors. 

It is eagerly desired to know the magnitude of mass itself, which would be given 
by the neutrinoless double beta decays. It is, however, impossible to know the 
neutrino masses at present. So, for the moment, it is important to pinpoint some 
of the mass-squared differences to the required accuracy by doing various experi- 
ments. One of the experiments, which is now running, is the K2K experiment [0] . 

In this paper, we investigate the oscillations of the beam, which is being 
measured at the K2K, in the four-neutrino model with mass scheme of the two 
nearly degenerate pairs separated by the order of leV for the three active and one 
sterile neutrinos [8- 15] by using the constraints on the mixing matrix derived from 
the solar neutrino deficit, atmospheric neutrino anomaly, Bugey reactor experi- 
ment, CHOOZ experiment, LSND experiments, CHORUS and NOMAD experi- 
ments and the other accelerator and reactor experiments. And, we compare the 
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oscillations with the ones in the three-neutrino model. 

Under the neutrino oscillation hypothesis ||16|| [|Hj], the flavor eigenstates are 
the mixtures of mass eigenstates with mass m 8 (i = 1, 2, 3, 4) in the four-neutrino 
model as follows: 

4 

Va = 52U 0fi v i , a = e,fx,r,s (1) 

i=l 

where i> e , and v T are the ordinary neutrinos and v s is the sterile one, and U 
is the unitary mixing matrix. The neutrino oscillation probability of u a — > v@ in 
vacuum is given in the usual manner by 

P (y a _> Vf) ) = £ Ke{U* ak U aj U* pj U pk ) sin 2 Ay+2 £ \m{Ul k U aj U; 3 Up k ) sin2A fej , 

k>j k>j 

(2) 

where A k j = Am| J L/(4£'), L being the distance from the neutrino source and E 
the energy of neutrino. 

The four neutrino masses should be devided into two pairs of close masses 
separated by a gap of about leV in order to accomodate with the solar and 
atmospheric neutrino deficits and the LSND experiments along with the other 
results from the accelerator and reactor experiments on the neutrino oscillation. 
There are the following two schemes for the mass pattern; (i) Am 2 olar = Am^ -C 
Am ltm = A «4 < A^lsnd = Am 2 2 , and (ii) Am s 2 olar = Am 2 3 < Am 2 tm ee 
Am^ Am 2 gND ee Am^, where Am^ ee m k — m 2 -. We will adopt the first 
scheme in the following analyses, and the second scheme can be attained only 
through the exchange of indices (1, 2) <-» (3, 4) in the following various expressions 
such as the oscillation probabilities. The constraints on the mixing matrix U are 
derived in the four-neutrino model in the following [|I~8| . 

(i) Solar neutrino deficit. Since A21 ~ 1 and all the other five A^'s are 
enormously larger than 1, the survival probability of u e is given from Eq.(|D by 



PsoUve^Ve) - 1 - 4|C/ el \ 2 \U e2 \ 2 sin 2 A 21 - 2\U e3 \ 2 (l - \U e3 \ 2 - \U t 



el | 
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-2\U e4 \ 2 (l-\U e4 \ 2 ), (3) 

where the unitarity of U is used. For the solar neutrino deficit, there are four 
different kinds of solutions as stated above, and a unique solution is not yet found, 
so that we will not use this deficit in order to obtain the constraints on U. 

(ii) Atmospheric neutrino anomaly. Since A 21 <C 1, A 43 ~ 1 and A 41 , A 42 , A 31 , 
A32 3> 1, the survival probability of is given by 



PatmK^^) ^ l-4|f/ M 3| 2 |^4rsm 2 A 43 -2(|[/ Ml | 2 + |f/ M2 | 2 )(l-|^ 1 | 2 -|^ 2 | 2 ). 

(4) 



By using the data from the Super-Kamiokande experiments, sin 2# atm > 0.82 for 



5x10 4 < Am 2 < 6x10 3 eV 2 , and expecting from this data that |Z7«i| 2 +|Z7«2| 2 <^ 



1, the following constraint is obtained, 

\V ie tf\V l $> 0.205. (5) 

(iii) The Bugey experiment [|19|1 (including Krasnoyarsk |[20|| , CDHS[21] and 
CCFR[^j experiments). By being typically represented by the Bugey reactor ex- 
periment with L/E = 3 — 20 [m/MeV], since A 24 <C 1, A 43 <C 1 and A 4i , A 42 , A31, 
A 32 ~ 1, the survival probability of z/ e is given by 

^Bugey^e - V e ) ~ 1 - 4( | U e3 \ 2 + | U e4 \ 2 ) ( 1 - | U e3 1 2 - | U eA \ 2 ) sin 2 A 32 . (6) 

If we use the data from the Bugey experiment conservatively, sin 2 26 ) Buge y < 0.1 
for 0.1 < Am 2 < 1 eV 2 , the following constraint is obtained: 

|t/ e3 | 2 + \U e 4 2 < 0.025. (7) 



The first long-baseline reactor experiment, that is, the CHOOZ experiment 
with L/E ~ 300 [km/GeV] gives a constraint of 4|?7 e3 | 2 |f/ e4 | 2 < 0.18 through their 
data of sin 2 2# C Hooz < 0.18 for 3 x 10~ 3 < Am 2 < 1.0 x 10~ 2 eV 2 . However, this 
constraint can be included in the constraint of Eq.(7) from the Bugey experiment. 



In the same way as the above, the LSND experiments^] with L/E = 0.5 — 1 
[m/MeV] brings the constraint of 

\U*^U e3 + U;M = 0.02 - 0.16 (8) 

from the data of sin 2 2^ LS nd = 1.5 x 10~ 3 -1.0 x 10" 1 for 0.2 < Am 2 < 2 eV 2 . And, 
CHORUS l24}| and NOMAD experiments searching for the — > z/ r oscillation 
with L/E = 0.02 - 0.03 [km/GeV] give the constraint of 

\u; 3 u T3 + u; 4 u T ,\< o.28 (9) 

from the data of sin 2#nomad < 0.3 for Am 2 < 2.2eV 2 . Therefore, among the 
abovementioned six typical phenomena and experiments, the useful constraints 
are of Eqs. (§, (g), @ and ®. 

In order to translate these four constraints to the ones for the mixing an- 
gles and phases, we adopt the most general parametrization of the mixing ma- 



trix for Majorana neutrinos, proposed by Barger, Dai, Whisnant and Young [|T5| . 
which includes six mixing angles and six phases. The expression of the matrix is 
too complicated to write it down here, so that we cite only the matrix elements 
which are useful for the following analyses; U e \ = C01C02C03, U e i = co2C^s* dQ1 , 

U e 3 = C 3S^ 02 , U ei = S^, U ^ = — Sd0 2 S d03«rfl3 + Co 2 Ci 3 S^ l2 , = C 03 S^ l3 , 

u r3 = 

C l3 S d02 S d03 s d23 C 02 s dl2 s dl3 s d23 + c 02Cl2C23; & n d U T 4 — Cq 3 Ci 3 S^ 23 , where 

= cos% and s dij = s^e" 5 ^ = sin%e 1<5lJ @, and 9 01 , 6 02 , 603, 612, #13, #23 are 



the six angles and 5 i, S 02 , S 03 , 5i 2 , 6i 3 , 623 are the six phases. Three of the six 
oscillation probability differences are independent so that only three of the six 
phases determine the oscillation probabilities, that is, the Dirac phases. 

By using this parametrization of U, the four constraints of Eqs. (5), (7), (8) 
and (9) are expressed by the angles and phases as follows: 

I - s 2So3Si 3 e- i(5o2 - 5o3+5l3) + c 02Ci 3 Si 2 e- i512 | 2 c 2 3S 2 3 > 0.205, (10) 
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44 + 4 < 0.025, (11) 

|C02S 2C03S1 2 C13 + cl 2 C 3So 3 Sl3e lSl 1 = 0.02 - 0.16, (12) 
c 02 C 12 s 12Cl3C23 — Co2So2So3-5i2Ci 3 S23e 1 ( <5l+<52 ) — Co 2 SQ 2 SQ3Ci 2 Si3C 2 3e lSl 

+ ci 3 si 3 s 23 (4 - cl 2 s\ 2 + 44)e~ i52 1 < 0.28, (13) 

where 5\ = 5 02 — <5o3 — ^12 + ^13 an d 8 2 = ^12 — ^13 + ^23- The constraint of Eq.(|l0|) 
reduces to 

s 2 12 c 2 13 s 2 13 > 0.205 (14) 

due to the smallness of S02 and S03, which is obtained from Eq.(ll). By using 
the constraints of Eqs.(ll) and (14) together with the nearly maximal mixing in 
the angle #23, which is derived from the large angle mixing in — > v T oscillation 
for the atmospheric neutrino anomaly, it proves that no constraints on the phases 
£1 and 62 are obtained from Eqs.(12) and (13). Equation (12), however, gives a 
constraint on the mixing angles. So, we obtain three constraints of Eqs.(ll), (12) 
and (14) on the mixing angles and no constraints on the two phases of 5i and 82 
in the four-neutrino model. The third Dirac phase does not give any significant 
effect to the leading parts of the oscillation probabilities, as can be seen from its 
no occurence in the constraints of Eqs.(10) — (13). 

By using the constraints of Eqs.(ll), (12) and (14), we calculate the oscillation 
probabilities of muon neutrino for the K2K experiment. We set the distance 
between the neutrino detector and the source L to be 250 km and the energy 
of neutrino E to be 1.4 GeV. A typical result is shown in Fig.l, which gives 
the survival probability Piy^ — > v^) and disappearance probabilities of Piy^ — > 
z/ r ),P(z/ At — > z/ e ), and P{y^ — > v s ) with respect to Am| 3 , which corresponds to 
Am 2 tm . We took the parameters as S01 = S23 = 1/^,^02 = S03 = 0.11, S12 = 
0.91, s 13 = 0.67, and 8\ = 8 2 = ir/2. Here and in the following, we take Am| 1 (= 
Am^ lar ) = 1.0 x 10" 6 eV 2 and Am^E Am^ SND ) = 0.3 eV 2 . The difference of 
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P( u fi ~~ > v ti) °f 0-7 from 1 around 0.0001 < AmJ 3 < 0.001 eV 2 comes from the 
^ m LSND contribution, as pointed out by Yasuda PBJ. It goes up to 1.0 — 0.96 in 
the same region of Amf 3 , as the angle S12 is increased towards 1.0. As seen in 
Fig.l, P{u^ — > Vy) does not vary in 0.0001 < Am^ < 0.001 eV 2 and decreases 
abruptly from 0.6 to 0.07 in the region of 0.002 < Amj 3 < 0.006 eV 2 . P{y^ u T ) 
and P(Ufj, — > u s ) are very small in 0.0001 < Am| 3 < 0.001 eV 2 and take a sizable 
magnitude of 0.3 — 0.6 in Am 2 3 = 0.004 — 0.006 eV 2 . In the computation we 
have not convoluted the probabilities with respect to the energy spread of the 
incident neutrinos. The feature of the computational results is as follows: (i) The 
dependence of the oscillation probabilities on the phases 5\ and 62 is very weak, 
(ii) The change of the probabilities between S02 = Sq^ = 0.11 and 0.05 is very 
small for Piu^ — > u^^Piu^ — > u T ) and Piu^ — > u s ), while the change is not so 
small for P{u^ — > u e ) since the mixing angles s 2 and s 3 affect significantly the 
v n ~^ v e oscillation, (iii) The dependence on the sign of the mixing angles Syi 
and S13 is that P{u^ — > u^) is same between the cases of su > 0, S13 > and 
S12 > 0, S13 < 0, and Piu^ — > u T ) and P(u^ — > z/ s ) interchange with each other 
between the two cases, and that all the probabilities do not change between the 
cases of Si 2 > 0, s 13 > and si 2 < 0, s% 3 < and also the same between the cases 
of S12 > 0, S13 < and s\2 < 0, S13 > 0. 

Next, we will compare these results with those in the three-neutrino model 
with three active neutrinos, where the mass pattern is taken as Am^ = Am 2 olar 
and Am 32 = Am 2 tm . The constraints on the mixing matrix are derived as follows. 
The solar neutrino deficit is not used again here, since there are so many (four) 
solutions. 

(i) Atmospheric neutrino anomaly. The survival probability of is given by 

PatmK -> Vy) ^ 1 "4|^ 3 | 2 (1 - | U & \ 2 ) Sin 2 A 32 . (15) 

Then, the constraint is replaced by | 1 2 ( 1 — |^3| 2 ) > 0.205, instead of Eq.(5). 

7 



This gives us the following constraint on (L^l, 

0.54 < < 0.84. (16) 

(ii) The CHOOZ experiment. The survival probability of v e is given by 

Pchooz(^ -> i7 e ) ~ 1 - 4|f/ e3 | 2 (l - |f/ e3 | 2 ) sin 2 A 32 . (17) 

The data stated in the description of the four-neutrino model gives the constraint 
of 4|[/ e3 | 2 (l — |t/ e 3| 2 ) < 0.18, leading to the following constraint on |[/ e3 |, 

\V*\ < 0.22. (18) 

In the three-neutrino model, the constraint from the Bugey experiment can be 
included in the constraint of Eq.(18) from the CHOOZ experiment. 

(iii) The CHORUS and NOMAD experiments. The transition probability of 
-> v T is given by 

P{v^v T ) ~4|f^ 3 | 2 |^ T3 | 2 sin 2 A 32 . (19) 

Their data of sin 2 26'nomad < 1 for Am 2 = (1.5 — 5) x 10~ 3 eV 2 gives a constraint 
of 

IMI^ral <0.5. (20) 

In order to translate these three constraints of Eqs.(16), (18) and (20) to the 
ones for the mixing angles and phases, we choose the parametrization used by 
the PDG[^7| with the additional two phases for Majorana neutrinos [|28[], which 



eventually includes three angles and three phases. We give here only the relevant 
elements; U e3 = si 3 e l(p_<w , = s 2 3Ci 3 e l(p ~ /3) , and U r3 = c 23 ci 3 , where c {j = cos% 
and Sij = sin 0^ and p, and f3 are the three phases. As in the four-neutrino model, 
the Majorana phases p and j3 do not enter into the oscillation probabilities. The 
three constraints are expressed with the three angles as 

0.54 < |s 23 ci 3 | < 0.84, (21) 



|s 13 | < 0.22, (22) 
c 2 3 |c 23 s 23 | < 0.5. (23) 

The constraint of Eq.(23) is always satisfied for any values of s 23 and s 13 . So, 
we have two constraints of Eqs.(21) and (22) remained. We again calculate the 
oscillation probabilities of the beam for the K2K experiment. A typical result 
is shown in Fig. 2 for Piy^ — > z/^),P(z/^ — > u T ), and Piy^ — > u e ) with respect to 
Am 2 2 (= Am 2 tm ), where we took s 12 = s 23 = 1/a/2, s i3 — 0.15, and = 7r/2. The 
probability -> is nearly 0.95 - 1 in 0.0001 < Am 2 2 < 0.001 eV 2 and 

decreases rapidly from 0.8 to 0.04 in the region of 0.002 < Am 2 ^ < 0.006 eV 2 . 
Piy^ — > v T ) takes a sizable magnitude of 0.2 — 0.6 around Am 2 2 = 0.002 — 0.004 
eV 2 and Piy^ — > f e ) is quite small due to the small mixing angle si 3 , which 
comes from the constraint by the CHOOZ experiment. The dependence of the 
probabilities on the phase is very weak. 

Fig. 3 shows the regions of Piy^ — > v^) allowed by the constraints on the mixing 
angles and phases discussed above in the three-neutrino model(dashed lines) and 
in the four- neutrino model(solid lines). The allowed region is the one sandwiched 
between the upper curve(maximum) and the lower one(minimum). It is very broad 
in the four-neutrino model, especially in the range of 0.0001 < Am 2 < 0.002 eV 2 
due to the range of 0.91 < s i2 < 1 derived from the constraint of Eq.(14), where 
the minimum curve of Piy^ — > v^) is given by s 12 = 0.91 and the maximum one 
is by s±2 = 1.0. On the contrast, the region of P{y^ — > v^) is considerably limited 
in the three-neutrino model, determined by the range of 0.55 < s 23 < 0.85 coming 
from the constraint of Eq.(21). The minimum curve is given by s 23 = l/y/2 and 
the maximum one is by s 23 = 0.55(or 0.85), where the maximum curve almost 
overlaps with the one in the four-neutrino model. The two models do not show 
any significant difference to Piy^ — > v^) in the range of 0.004 < Am 2 < 0.007 
eV 2 , while they give a large difference for Am 2 < 0.004 eV 2 . For example, if the 



9 



probability Piy^ — > v^) is measured to be 0.6 by the K2K experiment, the three- 
neutrino model will predict that the relevant scale of mass-squared difference is 
0.0030 < Am 2 < 0.0034 eV 2 , a relatively narrow range, while the four- neutrino 
model will predict that it is 0.0020 < Am 2 < 0.0034 eV 2 , a broader range than in 
the three-neutrino model. That is to say, the determination of the relevant scale 
of Am 2 by the K2K experiment is more flexible in the four-neutrino model than 
in the three-neutrino model. That is a big advantage at the time when Am 2 tm 
is not so precisely determined by the experiments and/or the observations, for 
example, of the atmospheric neutrino deficit. 

A way to predict more precisely the allowed region of Pivp, — > v^) for the K2K 
experiment in the four-neutrino model is to determine the angle s±2 to a good 
precision, which will be attained by the CHORUS/NOMAD-type experiment for 
the — > v T oscillation performed with a longer distance L ~ 30 km so that L/E 
becomes ~ 1, since Piy^ — > v T ) is predominantly controlled by the first term with 
the angle s±2 in the left-hand side of Eq.(13). 
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Figure captions 



Fig.l. The probabilities Piy^ — > v^) (solid line), Piy^ — > v T ) (dash-dotted 
line), Piyp — > z/ s ) (dashed line), and Piy^ — > z/ e ) (dotted line) versus Am| 3 
calculated in the four-neutrino model for the K2K experiment with L = 250 km 
and E = 1.4 GeV. The parameter values of the mixing angles and phases are 
soi = s 2 3 = l/V2,so2 = S03 = 0.11, S12 = 0.91, si 3 = 0.67, and Si = S 2 = n/2. 

Fig. 2. The probabilities Piy^ — > v^) (solid line), Piy^ — > i/ T ) (dash-dotted line), 
and -P(^ M — > f e ) (dotted line) versus Am^ calculated in the three-neutrino model 
for the K2K experiment. The parameter values of the mixing angles and phases 
are s i2 = S23 = I/a/2, S13 = 0.15, and = 7r/2. 

Fig. 3. The regions of Piy^ — > z^J for the K2K experiment allowed by the con- 
straints on the mixing angles and phases in the four-neutrino model (solid lines) 
and in the three-neutrino model (dashed lines). The allowed regions are in be- 
tween the upper and lower curves. The upper curve in the three-neutrino model 
almost overlaps with the upper curve in the four-neutrino model. The parameter 
values are s i = S23 = 1/a/2, S02 = 0.024, s 3 = 0.0, si 2 = 1.0, S13 = 0.54 and 
Si — S 2 — 7r/2 for the upper curve, and s i — $23 = 1/a/2, $02 = S03 = 0.11, Si 2 = 
0.91, S13 = 0.67 and Si = 0, S 2 = n/2 for the lower curve in the four-neutrino 
model. S12 = 1/a/2, S23 = 0.55, S13 = 0.20 and (p = ir/2 for the upper curve, and 
S12 = S23 = l/y/2, sis = 0.20 and = 7r/2 for the lower curve in the three-neutrino 
model. 



14 



0.0 001 0.0 01 0.01 

Am 2 (eV 2 ) 

32 

Fig .2 



